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Abstract 

The purpose of our study is to understand the mathematical origin in real 
space of modulated and damped sinusoidal peaks observed in cosmic mi- 
crowave background radiation anisotropies. We use the theory of the Fourier 
transform to connect localized features of the two-point correlation function 
in real space to oscillations in the power spectrum. We also illustrate ana- 
lytically and by means of Monte Carlo simulations the angular correlation 
function for distributions of filled disks with fixed or variable radii capable 
of generating oscillations in the power spectrum. While the power spectrum 
shows repeated information in the form of multiple peaks and oscillations, 
the angular correlation function offers a more compact presentation that con- 
denses all the information of the multiple peaks into a localized real space 
feature. We have seen that oscillations in the power spectrum arise when 
there is a discontinuity in a given derivative of the angular correlation func- 
tion at a given angular distance. These kinds of discontinuities do not need 
to be abrupt in an infinitesimal range of angular distances but may also be 
smooth, and can be generated by simply distributing excesses of antenna 
temperature in filled disks of fixed or variable radii on the sky, provided 
that there is a non-null minimum radius and/or the maximum radius is con- 
strained. 

Keywords: cosmic microwave background, Methods: statistical 



Preprint submitted to Physica A 



September 21, 2012 



1. Introduction 



We present a series of two papers with the motivation of understand- 
ing the CMBR ("Cosmic Microwave Background Radiation") anisotropics [l| 
without direct comparison with a physical model, in order to see whether 
other conditions different from the standard model of cosmology may gener- 
ate something similar of what is observed. In this first part, we will focus 
qualitatively on the mathematical elements necessary to produce peaks in 
the power spectrum of the CMBR, and in the second part [2j we will discuss 
its physical generic interpretation and evaluate quantitatively the number of 
free parameters a generic theory needs in order to fit it. This line of research 
is at present almost totally abandoned, the available literature merely dis- 
cussing how the peaks can be generated in terms of the standard cosmology. 
We could consider our theme as an inversion problem: rather than deriving 
the predictions from a model to fit the data, we wish to ascertain from the 
data some general characteristics of a generic model. Our goal is to open 
analysis to wider mathematico-physical scenarios, in order to see whether 
something different from standard cosmology interpretation can reproduce 
it. 

The analysis of CMB maps is usually performed in terms of spherical har- 
monic decomposition and by computing the angular power spectrum of CMB 
anisotropies pi, 0, |i[ Q . The standard cosmological model provides a phys- 
ical model able to fit accurately the power spectrum of CMB anisotropies. 
It is well known that this exhibits a series of peaks, or an almost periodic 
oscillation, as predicted by Peebles & Yu Q. 

In the ideal case of full-sky coverage, the two-point angular correlation 
function provides a complementary means of analyzing CMB observations 
instead of the angular power spectrum and, in principle, contains the same 
information as the angular power spectrum. However, as we will discuss in 
what follows, the angular correlation function allows an easier understanding 
of the anisotropy structures and it may serve as a complementary means of 
spherical harmonic analysis. Few authors have considered the direct deter- 
mination of the anisotropies correlation function directly in angular space 
jsl, 0, [lo|, 11, 12, 13, 14], (see the review in ref. 15]), and many of them 



have focused on the large scale behavior of the angular correlation function, 
where the effect of cosmic variance is certainly important. In particular, the 
two-point correlation function almost vanishes on scales greater than about 
60 degrees, contrary to what the standard CDM theory predicts, and is in 
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agreement with the same finding obtained from COBE data about a decade 
earlier jgl, It was pointed out[15| that the striking feature of the two-point 
angular correlation function is the lack of large angle correlations, unexpected 
in inflationary models. On the other hand it is also interesting to consider, 
both theoretically and observationally, the behavior of the correlation func- 
tion on small angular scales 16lll7|. Indeed, the whole structure of the peaks 
in the angular power spectrum corresponds in direct angular space to a lo- 
calized feature of the correlation function which is analogous for the CMB 
field to the baryon acoustic oscillation (BAO) scale characterizing the matter 
correlation function 18, [l9[ (see criticism on its detection in ref. |io|)- The 



direct identification of such a scale would thus represent a complementary 
test of the standard model. 

In this paper, we illustrate some simple mathematical properties of the 
anisotropy structures associated with both the oscillating peaks in the an- 
gular power spectrum and with the characteristic scale in the anisotropy 
angular distribution. In particular, we provide a simple explanation of the 
peak sequence, showing that it is simply associated with the shape of the 
structures of the anisotropies. While this discussion is useful to illustrate 
the statistical meaning of the properties of the A-CDM angular correlation 
function, it also provides a simple general model to generate such type of 
correlations. 

Our paper is structured as follows: In §0 we discuss the mathematical 
origin of the oscillating peaks in the power spectrum, and of the correspond- 
ing scales in the angular distribution of anisotropies. In £j3l we illustrate 
the situation by introducing a few simple toy models able to show similar 
correlation properties as those predicted by standard cosmological models. 
Finally, in §3]we summarize the results. 



2. The origin of oscillations in the power spectrum and its corre- 
spondence in the correlation function at small angles 

2.1. Power spectrum and self- correlation 

The fluctuation field on the sky, ^r-(n), where n is the general point on the 
unit sphere (or direction in 3— dim space), can be decomposed into spherical 
harmonics in the sphere: 

^( n ) = Yl a ^^ m (n) . (1) 
1=0 m =-i 
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Assuming statistical isotropy, the angular power spectrum Ci is defined as 



Ce = (|^ m | 2 ) 



2£ + 



\ a £r, 



(2) 



m=—t 



whereas the angular two-point correlation function is defined as [ 21] . [22| (§6.6.1)]: 

C(9) = ^(n)^(n')^ (3) 

where 9 is the separation angle between the two directions n and n'. The two- 
point correlation function can be expanded in terms of Legendre polynomials 



By inverting Eq. (J4]) we find [151 



c(e) = — ^(2e + i)c t p t [co8 



(4) 



C £ = 2n / C(9)P e [cos(9)\ sin(0)d0. 
'o 



(5) 



Alternatively, always assuming statistical isotropy, we can express the 
power spectrum in terms of the conjugate frequency k in the Fourier trans- 
form: 



P(k) = C(k) = FT[C{6)} = / d9C(9)e 



-ik9 



(6) 



where the vectors 9 and k are two-dimensional vectors whose moduli are 9 
and k respectively, and 



C(9) =FT~ 1 [P(k)] . (7) 

For small enough angles < 1 radians, P e [cos(6)] « J [(£ + 1/2)9} ^(Eq. 
46.39), where Jo is the Bessel function of order zero, and P(k) ~ 2tx J Q n C{9) 
Jo(k9) sin 6 ] d9 [23j(Eq. 46.43), so the spherical harmonic decomposition co- 
incides with the two-dimensional Fourier transform with a frequency k = 
[I + 1/2). From now on, we will not distinguish between Ce and P(k). 
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2.2. Relationship between oscillations in the power spectrum and abrupt tran- 
sitions in the self- correlation function 

In order to clarify the relation between the sequence of oscillating peaks in 
the CMB angular power spectrum and the angular space two-point correla- 
tion function, let us introduce a simple model. We first distribute points cen- 
ters uniformly in space and then distribute continuous mass profiles around 
these points, which hereafter we denominate as centers. Since we are in- 
terested in the angular correlation function C{6) at 9 <C 1 radians, we can 
approximate the angular space as Euclidean and use Fourier transform analy- 
sis to illustrate mathematically the link between sinusoidal oscillations in the 
power spectrum and strongly localized features in the two-point correlation 
function. 

In general in a d— dimensional Euclidean space the model can be formu- 
lated as follows. Let us start with a stochastic spatial distribution of centers 
with microscopic number density 



where the sum runs over all center positions Xj. Let us assume that this 
distribution is statistically homogeneous in space, which implies: (i) the 
average number density is a positive constant n$ > 0, and (ii) the two-point 
correlation function (TPCF) depends only on the separation vector between 
these two points. In this case, the TPCF can be defined as [22[ 




where (•) represents the ensemble average, or the infinite volume average in 
the case of ergodicity. We assume that h(x) (the covariance or off-diagonal 
correlation function) is a sufficiently regular function with no discontinuity 
in any of its derivatives for x > and decays sufficiently rapidly at large 
separations x (e.g., one can suppose an analytic function in all space x as an 
exponential or power law decay). 

Due to the assumed statistical and spatial homogeneity, the power spec- 
trum S(k) of the stochastic distribution of centers satisfies 




(8) 







(10) 
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where 



Sh(k;V)= [ d d x[n(x) - n ]e- ikx (11) 
Jv 

h(k) = FT[h(x)] = [ d d x /i(x)e- ik ' x . (12) 

Note that, since h(h) vanishes at large k, S(k) in the same limit converge 
to 1. As shown more rigorously below, since h(x) is thought to have all 
continuous derivatives at x > 0, S(k) is expected not to present damped 
sinusoidal oscillations. 

Let us now replace each center at Xj with a continuous density profile 
described by the function /(x— Xj) > 0, such that /(x — > oo) — > sufficiently 
fast so that J Rd d d x f(x) is equal to a finite positive constant which we can 
take to be 1. The field can now be written as 

T(x) = ^)/(x-x i ). (13) 

i 

The mean value is clearly (T(x)) = n$ j d d x /(x). 
Since we can write 

T(x) = J d d yf(x-y)n(y) } (14) 

it is simple to show that the power spectrum -P(k) of the field T(x) is simply 
given by 

P(k) = |/(k)| 2 S(k), (15) 

where /(k) = FT[/(x)]. Note that, as /(x) is assumed to be integrable, 
/(k) vanishes at large k and therefore also -P(k) vanishes at large fc. At the 
same time, the TPCF of the field T(x) is a continuous function, as it has to 
be for any proper continuous field |22|. 

Since we have assumed that S(k) is not characterized by oscillating peaks, 
if we want -P(k) to show them, we need these to be associated with /(k). 
Thus, the problem is reduced to finding which features of /(x) generate 
oscillating peaks in -P(k). We show below that if /(x) presents a finite 
discontinuity in the function or in one of its derivatives at a position x ^ 0, 
then -P(k) will show slowly damped sinusoidal oscillations and wave-vector 
peaks (in Fourier space) x . 
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A simple example is given by the following function /(x) in d dimensions, 
usually called spherical box function 



/(*) 



l/\\S d (R)\\ foixeS d {R) 

otherwise 




where Sd{R) is the d— dimensional sphere of radius R centered at the origin 
and ||/Sd(i2)|| its volume. In this case, /(k) is characterized by slowly damped 
sinusoidal oscillations: 

sin(fci?) 
~kR ' 

/(*) = < TZaMkR), d = 2 } , (16) 

o sin(fc_R)— kRcos(kR) 

where J% is the Bessel function of the first kind. If, for example, we have a 
Poissonian distribution of centers, i.e., one with S(k) = 1, we have simply 
-P(k) = \f(k)\ 2 which shows slowly damped sinusoidally oscillating peaks. 
A generalization of this development for a variable radius R is given in the 
following subsection. 

2. 3. Mathematical description of the origin of oscillations in the varying ra- 
dius case 

We can generalize the mathematical description of §2.2l for variable radius 
balls by considering the following field 

T(x) = ^/(x-x i ;i2 i ), (17) 

i 

where {xj} are the positions of the centers and 

/(x; R i )=A(R i )g(x)9(R i -x), 

where the radius Ri changes from center to center, g(x) is an arbitrary uni- 
versal non-negative function and A(Ri) is the probability density function. 

In order to find the power spectrum P(k) we should calculate the following 
double average 



\ST(k,V)[- 
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where 5T(k,V) = FTy[T(x) — T ], (•) is the ensemble average over the 
positions of centers and (•) is the average over the distributions of radii. The 
two averages clearly commute. The first step is to write 



T(k,V) 



d d xT(x.)e 



We can then write 



\f (k, v)\ 2 = />, Ri)ffc R 3 )e- l ^-^ ] • 



It is simple to show that the diagonal part of this double sum is exactly 
balanced by the subtraction of the mean density from T(x) in the actual 
definition of the PS, for this reason only the non-diagonal part contributes 
to the final PS. We have therefore to calculate the following average 



where the prime in the double sum means the condition i ^ j, and 



Thus the problem is reduced to calculating /(k; R). 

In other words, the power spectrum is the superposition of the power 
spectra for each value of Ri. 

2. 4- General relation between singular points in a function and damped os- 
cillations in its FT 

Let us now give the general argument relating the discontinuity of the 
function /(x), or of one of its derivatives at some point xo 7^ to the pres- 
ence of slowly (i.e., power law) damped oscillations (and peaks) in P(k). 




(18) 




By taking the second average (■) of Eq. ( IIS]) , finally we get 



P(k) = f(k;R) S(k). 
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With the aim of simplicity let us develop the argument in d = 1 dimen- 
sions. Let us assume that the function f(x) is analytic everywhere with the 
exception of the point xq where its n th derivative f^ n \x) presents a finite 
discontinuity: /^(xq ) — /^(xq ) = A. Moreover, let us assume that f(x) 
vanishes sufficiently rapidly at large \x\. Let us study the features of the FT 
of f(x) induced by such discontinuities. The FT is defined by 

/+oo 
dxe~ lkx f(x). (19) 
-oo 

By integrating it by parts n + 1 times, we get 

(20) 



where 

g(x) = f {n+1 \x) - AS(x - x ) 

is an analytic function everywhere but in xq where at most it can have a 
finite discontinuity. In this last case this discontinuity will lead to a further, 
but subdominant (i.e., decreasing faster with k) contribution to the slowly 
damped oscillations given by the last term in Eq. f[2"Uj) . It is important to 
note that the power of decay of damped oscillations is directly related to the 
order of the discontinuous derivative. 

This argument can be extended to higher dimensions. For a statistically 
isotropic distribution in d dimensions in which both /(x) and £(x) depend 
only on x = |x|, and consequently /(k), S(k) and -P(k) depend only on 
k — |k|. In this case, when /(x) = f(x) shows a discontinuity in only one 
of the derivatives in x, the calculation can be reduced to a one-dimensional 
one very similar to the one-dimensional case treated above. 

The terms exp(— ikx ) in Eq. ( 120]) give the oscillations. Thus oscillat- 
ing peaks in the power spectrum occur when the correlation function has an 

9 



{ik) n+1 



, —, lim 

{lk) n+l e,6'^0 



oo 
xo-e 



dxe~ ikx f {n+1) {x) 
dxe- ikx f {n+1 \x) + 



dx e~ ikx f n+1 \x) + A \ dxe~ ikx 5{x 



x +e' 



xo+e' 



xo—e 



1 



(ik) n+l J 



dxe x g(x) + 



j^ e ~ikx 
(ik) n+1 



abrupt change of functionality at some 9 > 0; that is, when there is a char- 
acteristic angular distance which separates two regimes in the self-correlation 
C(9); mathematically, this means that some derivative is not continuous at 
this point. For our case, d = 2, this happens in any TPCF such that 

C(0) = I 9 A 6 } lz H . (21) 



92(9), 9>9 







with^^^W)- 

The exact structure of the oscillations depends both on the point with 

some non-continuous derivative, 9q, and also on the slope before and after 

that point: gi(0Q ) and q 2 {9q). In general, it may also happen that there is 

more than one point with a point of this kind. In the case of CMBR analysis, 

one point with some non-continuous derivative is enough to explain all the 



3. Examples and toy models 

In order to make explicit the discussion in Sect. |2] let us consider the 
following two TPCF C x (9) and C 2 {9): 

C x {9) = A 1 , 1 exp(-6/e 1>1 ) + A 1)2 exp(-0/0 1>2 ) , 

where we have chosen the following numerical values: A\ t \ = 9744, A\ >2 — 
3000,6*1 1 = 0.45,6*1,2 = 13.0 (angles are in degrees), and 

n(ff\-l ^2,iexp(-^ 2 , 1 ), 6 < 9, 
° 2W \ A 2)2 exp(-^ 2)2 ), 9>9* ' 

where A 2>1 = 12000, A 2:2 = 3600, 9 2>1 = 0.79,9 1>2 = 11.45,0, = 1.03. 

In Fig. [TJ we plot these two correlation functions and their corresponding 
power spectra, derived through Eq. (jSJ). 

While both correlation functions present a relatively fast change in be- 
havior at 9 1°, only the power spectrum of C 2 (9) presents an oscillating 
series of peaks. 

This simple example clearly shows that only the Fourier transform of a 
function which has a discontinuous derivative (i.e., C 2 {9)) presents a series 
of oscillating peaks. 

Let us now consider two illustrative toy models, constructed following the 
recipe discussed in Sect. 12.21 
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9 (deg.) 




Figure 1: Upper panel: Test correlation functions C\(9) and 02(6). Bottom panel: The 
corresponding power spectra FT[Ci] and FT[C2\- 
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Figure 2: Representation of the geometrical relationships between two points i, j at two 
disks. 

3.1. Toy model 1: disks with same radius 

We construct a statistically isotropic temperature field T(n) as follows. 
We distribute disks of angular radius R such that for separation angles from 
the center of each disc 9 < R the temperature field is T(n) = f(9 r ), 9 r being 
the distance from the center of the disk and the value of the field outside the 
disks being T = 0. There are N c disks in a given angular region; their centers 
are distributed as a spatially homogeneous stochastic point process with an 
angular two-point correlation function u(9). When oo(9) = for 9 > the 
centers of the disks are distributed as a homogeneous spatial Poisson process. 
Given that we are interested in low values of 9 <C 1 (radians) , we can use the 
approximation sin(0) ~ 9. 

The self-correlation of the field T corresponds to the integration over all 
possible pairs of points respectively within two disks (see Fig. [2]). The 
self correlation is the average of the product of the value of the temperature 
in the first disk, f{9j), multiplied by the value of the temperature in the 
second disk, f(9j). The cases in which some point i or j is not within some 
disk give a null contribution to the integral. Therefore, our problem reduces 
to constraining the limits of the integral of the average correlation within 
the cases with non-null contribution. The problem is not trivial but, after a 
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patient analysis, one can see that the solution is given by: 



C{9) 



An 2ix9 




d(f), 



(22) 



R 



X 



de t 9j(e t )[i s (e,e t ) + 1.(9,9,)}, 



where I s d9 and I d9 stand for the integrals of the value of T over all pixels 
Xj such that 9 < \xi — Xj\ < 9 + d9 respectively in the areas of the same 
circle of Xi (J s ) and in the areas of other disks (I Q ). They do not depend on 
(pi because both / and u depend only on the angular distance. 

I s is an integral of the second areas in the case that i and j are in the 
same disk, integrating over all possible values of 9j, (fij that keep the distance 
between the points i and j between 9 and 9 + d9. It is: 



The evaluation of 9j is a function of 9, 9i, <pj, such that \xi — Xj\ = 9: 



If, within the constraint, there are two possible values of 9j, both will con- 
tribute to the integral in Eq. ( )23|) . while if there are no values of 9j no value 
will be included in the integral. This is a special case of Eq. ( I2"gj) for 9 Q = 0. 

For the evaluation of the second integral involving I Q , we have to do 
something similar but moving the center of the disk j among all possibilities. 
The center of the "other" disk is at angular distance 9 Q from the first one, 
and azimuth <ft . Thus, 




(23) 



9j{9, Si, <pj) = 9, cos(0 j ) ± VA 
A = 91 cos(^) 2 - 91 + 9 2 
where A > 0; < 9 3 < R . 



(24) 




(25) 
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2tt 

d6 9 P(9 ) 



X 



d9 j {9,9 i , ( j)i,0oAo) 



06 



d9 



where P{9 ) is the probability per unit area of finding a circle at distance 9 Q 
from the center of the first circle. 9j is a function 9, 9i, <ftj, 9q, <p standing 
again for the value of 9j, which follows \xi — Xj\ = 9 but taking into account 
that 9j and <f>j are coordinates with respect to a second circle whose center is 
at coordinates 9 Q , O with respect to the first one. In addition, by definition 
we have: 

N 

P(9 ) = ^[l+u(9 )}. (26) 

For the calculation of 9j, by considering the two-dimensional planar ap- 
proximation (instead of spherical geometry) illustrated in Fig. [2] we find 

9 = \xi — Xj\ , (27) 

X{ 9ji , 

Xj = 9 Q (cos 4> Q i + sin (f> j) + ^(cos fyi + sin <pjj) . 

Hence, 

9j(9, 9i, fy, 9 Q , 4> ) = -9 a cos(0 o - <f>j) + 9i cos^j) (28) 

±Va 

A = [B cos(0 o - <P 3 ) - Bi cos^)] 2 ~ e l ~ ^ + Q 2 

+29A cos(0 o ) 

where A > 0; < 9 3 < R . 

As in the previous case, if there are two possible values of 9j, both of them 
contributing to the integral in Eq. ( 125]) . while if there are no values of 9j 
within the constraint, no value will be included in the integral. 

The results of this toy model for some test functions are given in Fig. [3] 
respectively for: 

a) f{9) = 1 (homogeneous distribution within the disk), uj{9) = (Pois- 
sonian distribution of disks), R = 1°, iV c = 10 3 ; 
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b) f(9) = 1, and 



u>{9) 







i, e < 2R 

6>2R 



(29) 



i.e., a Poissonian distribution, with the added restriction that disks do not 
intersect, R = 1°, N c = 10 3 ; 

c) f(6) = 1, u(6) = 2 x exp(-6/R) -1,R= 1°, N c = 10 3 ; 

d) f{8) = exp(-6/R), u{8) = 2exp(-0/R) - 1, R = 1°, N c = 10 3 . 

We have also carried out Monte Carlo simulations: We have randomly 
distributed N c centers, with or without the constraint that minimal distance 
between two nearest neighboring centers is larger than 2R, where R is the 
disk's size. Inside each disk of size R we randomly distribute N p points. 
When disks are in a non-overlapping condition the minimal distance between 
the N c centers has to be > 2R. We have chosen R = 1° and we distribute 
N c = 80 disks in a square of side L = 1 rad., in addition we fixed iV c = 1256; 
the correlation functions are computed as an average over 50 realizations of 
the distribution. The results of the analytical calculation of cases la) and 
lb) were compared with these Monte Carlo simulations (with the appropriate 
normalization to get the same value of f(8 < 1°) = 1). The results are plotted 
in Fig. |3J numerical simulations nicely agree with analytical calculations. 

3.2. Toy model 2: disks with varying radii 

With regard to disks of different size, we consider only the case in which 
oj{9) = 0, i.e., a Poisson distribution of disk centers. In the more general 
case of to (8) ^ 0, the shape of C{6) would be changed, as we have seen in 
the previous subsection. However the feature corresponding to the abrupt 
transition [i.e., a discontinuous first derivative of C(8)], which is the scale of 
interest for the generation of the oscillating sequence of peaks in the angular 
power spectrum, does not depend on the particular shape of uj{8). We are 
going to show that the angular TPCF of the set of disks of different size is 
again characterized by a discontinuous derivative. 

Let us assume that we have disks of different angular sizes R, and that, 
conditioned by their size, the contribution to the angular correlation function 
of disks with radius between R and R + dR follows a dependence 




(30) 
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analytical 
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analytical 
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- Case lb 


simulation 





) (deg.) 



Figure 3: Upper panel: Self-correlation C(9) for "toy-model 1" cases: a) f(9) = 1, Pois- 
sonian distribution of disks; b) as in case a but with the added restriction that disks do 
not intersect); c) as in case a with two-point angular correlation function of the centers 
of disks uj(9) = 2 x exp(—9/R) — 1; d) as in case c but with f{9) — exp(—9/R). All the 
cases with a number of disks N c = 10 3 of radius R = 1°. Bottom panel: Comparison of 
cases la) and lb) derived analytically with the results of Monte Carlo simulations; with 
the error bars, we give the r.m.s. of the simulation lb); for case la) is similar. 



16 



where A(R) is the amplitude of the correlation of disks with radius R, and h is 
a smooth continuous function with all continuous derivatives such that h{2) = 
0; such a distribution of h is given by the contribution to the correlation of 
the disks with the same size (i.e., the toy model discussed in the previous 
section); we subtract the average field [thus, the constant C{6 > 2R) may 
be set as zero conveniently defining a new T* = T — constant}. Assuming a 
uniform distribution of radii between -R m i n and -R ma x, given that there is no 
correlation between the positions of the disks, the angular TPCF is just the 
sum of the correlation function between its minimum and maximum radii. 
That is, 



C(6) = / C(9,R)dR = 

0fi!£ZdxA(e/x)h(x)/a*, 

9 fe/R max dxA (8/x)Hx)/x 2 , 
0, 



(31) 



@ — ^Rmin 
^Rmin <L 2R n 
@ ^ 2Rm.n.T. 



If there are minimum and maximum radii there is clearly a change of de- 
pendence with 9 for three regimes. If, instead, there is only a minimum 
{Rmin > 0) or a maximum radius, there will be two regimes. Finally, if 
Rmin = 0, Rmax = oo, then there will be only one regime and no abrupt 
transition. 

Let us consider an example: a) A = 1, h(x) = (1 — x/2), for which we get 



{Rmax R 
Rmax 

0, 



C{6) = 

2^{Rmax I Rmin)® i @ — ^Rr, 



(32) 



1 -±^e + Un{e/R max ) 1 



^ ^RmciT 



Fig. m displays Eq. flMJ for R min = 1° and R max = 2°. 

To each disk radius projected on the sky we could associate the physical 
distance R = L/r, where r is the distance and L is the physical size in an 
Euclidean space. As it will be discussed in jsjj], this may be representative 
of spheres placed at different distances. For A{R) = (^) and for instance 
h{x) = (1 — x/2), we find 



C{9) 



(33) 
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Figure 4: Upper panel: Two-point correlation C(6) for the disks of varying radius, when 
no correlation among disk centers is considered a), b). Bottom panel: The corresponding 
power spectra. 
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\ r max -r mzn LP — r 2 R < 

I L 4L 2 \max 'min)°\i u — u l 

[-^ + 1 + ^0], 6 1 <6<6: 



0>02 



9 X = 2L/r max ; 6 2 = 2L/r 



Fig. HJdisplays Eq. (1331) for a case b) v4 = 0.02, r min / L = 3, r max /L = 50. 

As can be observed, the correlation function is smoother. However, the 
power spectra (Fig. H|) still show a sequence of oscillating peaks. Indeed, 
another way to understand the fact that oscillations are still present when 
the disk size is not fixed is obtained recalling that the power spectrum is the 
superposition of the power spectra for each radius (see §2.3p : the superposi- 
tion of oscillating power spectra gives still a global oscillation, provided they 
do not cancel each other out (which is the case when there is no minimum 
and/or maximum radius). 

The examples of this subsection are some particular ways of convolution. 
In general, any kind of convolution or smoothing in the abruptness of the 
disc shape will also produce oscillations, though with damped amplitude of 
the oscillation. 

4. Conclusions 

Although the power spectrum analysis is useful for obtaining some kind of 
information on CMBR anisotropy data, the angular TPCF (Two-point corre- 
lation function) analysis gives a more straightforward physical representation 
of that distribution. Moreover, while the power spectrum shows repeated in- 
formation in the form of its multiple peaks and oscillations, the TPCF, its 
Fourier transform, offers a compacter presentation which condenses all the 
information of the multiple peaks in a localized real space feature. Precisely 
because of that, and because there is a dearth of literature analyzing the 
properties of the TPCF, we have concentrated here on its analysis. 

In this paper we have shown how oscillations in the power spectrum arise 
when there are some point singularities in the angular TPCF in the form of 
some discontinuity point of its n-th derivative at some angular distance. In 
particular, we have clarified the analytical link between the frequency and 
the power law order decay of damped oscillations in the power spectrum and 
the kind of singularity in the TPCF. We have then presented toy models 
able to generate and illustrate this phenomenology: simply by placing some 
distribution in the sky of filled disks of fixed or variable radii with an excess of 
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antenna temperature, provided there is a minimum> and/or a maximum 
in the range of radii. Disks of different radii may stand for projection of 
spheres with variable radii or spheres with variable/fixed radii and variable 
distance from the observer. 

Discussions on the physical interpretation of these mathematical prop- 
erties for the case of CMB anisotropies and comparison with real CMBR 
data will be given in j^j. This paper, now in preparation, will discuss these 
mathematical properties in terms of matter distribution in the fluid which is 
generating the radiation. The angular correlation function of the anisotropies 
from CMBR data will also be calculated, in an effort to derive the minimum 
number of parameters of a generic function to fit it. The paper 0] will argue 
that a power spectrum with oscillations is a rather normal characteristic ex- 
pected from any fluid with clouds of overdensities that emit/absorb radiation 
or interact gravitationally with the photons, and with a finite range of sizes 
and distances for those clouds. It will also show that the angular correlation 
function can be fitted by a generic function with a total of ~ 6 free parame- 
ters. The standard cosmological interpretation of "acoustic" peaks is just a 
particular case; peaks in the power spectrum might be generated in scenarios 
that have nothing to do with oscillations due to gravitational compression in 
a fluid. Nonetheless, the standard model with six free parameters produces 
a better fit than the generic fit with the same number of free parameters; it 
fits third and higher order peaks whereas the generic fit reproduces only the 
first two peaks. 

Apart from the analysis of CMBR peaks, the concepts analyzed in this 
paper are general and applicable to any kind of power spectrum with oscilla- 
tions/peaks. The BAO oscillations in the large-scale structure 19|] might be, 
for instance, another field where these results can be applied. 
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